/*	This program produces Figure D.7 */

***** Set directories 
local dir_clean 	"~/Dropbox/Retirement gaming/clean"
local dir_output 	"~/Dropbox/Retirement gaming/output/dataverse"


use "`dir_clean'/ech_clean.dta", clear

local listcontrols = " i.t i.fsize i.ciiu2  i.year#i.fsize i.year#i.ciiu2   " 


** SELF EMPLOYED NET EARNINGS
preserve
keep if self_empl==1
*Normalize hours
sum Wnet if age_centered==-1
local meanw=r(mean)
g Wnetz=Wnet/`meanw'
* Age dummies
tab age_centered, gen(age_centered_dums)
drop age_centered_dums5 
* Regression
reghdfe Wnetz age_centered_dums* [aw=weight], absorb(`listcontrols') vce(robust)
* Sample size
global n=e(N)
* Coefficients and sd
for any beta sd: gen X = .
forvalues d = 1(1)13 {
	capture qui replace beta = _b[age_centered_dums`d'] if age_centered==`d'-6
	capture qui replace sd   = _se[age_centered_dums`d'] if age_centered==`d'-6
}
replace beta=0 if age_centered==-1
replace sd=0 if age_centered==-1
* Plot the results
collapse beta* sd*, by(age_centered)
gen sd_top = beta + 1.96*sd	
gen sd_bot = beta - 1.96*sd
gen sd_top10 = beta + 1.645*sd	
gen sd_bot10 = beta - 1.645*sd
sum sd*
twoway rcap sd_top sd_bot age_centered, ///
	cmissing(n) lwidth(thin) lcolor(navy) lpattern(solid) xsc(r(-5 7)) xlabel(-5(1)7)  xline(-1) xline(3.5, lp(dash) ) ysc(r(-0.15 .15)) ylabel(-0.15(0.05).15,grid) ///
	|| rcap sd_top10 sd_bot10 age_centered, cmissing(n) lwidth(thin) lcolor(navy) lpattern(solid) xsc(r(-5 7)) ///
	|| scatter beta age_centered, mcolor(navy) scheme(s1color) lcolor(navy) lpattern(solid) yline(0, lcolor(gs2)) ///
	title("Earnings, relative to age-49 average")  note("N=${n}") ///
	legend(off) ytitle("Estimated Coefficients") xtitle("Years Relative to Age 50")   
graph export "`dir_output'/figureD7a.png",  replace 	
restore

		
		
** EMPLOYED FIRM<10 NET EARNINGS
preserve
keep if  empl==1 & small==1
*Normalize wages
sum Wnet if age_centered==-1
local meanw=r(mean)
g Wnetz=Wnet/`meanw'
* Age dummies
tab age_centered, gen(age_centered_dums)
drop age_centered_dums5 
* Regression
reghdfe Wnetz age_centered_dums*  [aw=weight], absorb(`listcontrols') vce(robust)
* Sample size
global n=e(N)
* Coefficients and sd
for any beta sd: gen X = .
forvalues d = 1(1)13 {
	capture qui replace beta = _b[age_centered_dums`d'] if age_centered==`d'-6
	capture qui replace sd   = _se[age_centered_dums`d'] if age_centered==`d'-6
}
replace beta=0 if age_centered==-1
replace sd=0 if age_centered==-1
* Plot the results
collapse beta* sd*, by(age_centered)
gen sd_top = beta + 1.96*sd	
gen sd_bot = beta - 1.96*sd
gen sd_top10 = beta + 1.645*sd	
gen sd_bot10 = beta - 1.645*sd
sum sd*
twoway rcap sd_top sd_bot age_centered, ///
	cmissing(n) lwidth(thin) lcolor(navy) lpattern(solid) xsc(r(-5 7)) xlabel(-5(1)7)  xline(-1) xline(3.5, lp(dash) ) ysc(r(-0.15 .15)) ylabel(-0.15(0.05).15,grid)  ///
	|| rcap sd_top10 sd_bot10 age_centered, cmissing(n) lwidth(thin) lcolor(navy) lpattern(solid) xsc(r(-5 7)) ///
	|| scatter beta age_centered, mcolor(navy) scheme(s1color) lcolor(navy) lpattern(solid) yline(0, lcolor(gs2)) ///
	title("Earnings, relative to age-49 average") note("N=${n}") ///
	legend(off) ytitle("Estimated Coefficients") xtitle("Years Relative to Age 50")   
graph export "`dir_output'/figureD7b.png",  replace 	
restore



** EMPLOYED FIRM>=10 NET EARNINGS 
preserve
keep if  empl==1 & small==0 
*Normalize wages
sum Wnet if age_centered==-1
local meanw=r(mean)
g Wnetz=Wnet/`meanw'
* Age dummies
tab age_centered, gen(age_centered_dums)
drop age_centered_dums5 
* Regression
reghdfe Wnetz age_centered_dums* [aw=weight], absorb(`listcontrols') vce(robust)
* Sample size
global n=e(N)
* Coefficients and sd
for any beta sd: gen X = .
forvalues d = 1(1)13 {
	capture qui replace beta = _b[age_centered_dums`d'] if age_centered==`d'-6
	capture qui replace sd   = _se[age_centered_dums`d'] if age_centered==`d'-6
}
replace beta=0 if age_centered==-1
replace sd=0 if age_centered==-1
* Plot the results
collapse beta* sd*, by(age_centered)
gen sd_top = beta + 1.96*sd	
gen sd_bot = beta - 1.96*sd
gen sd_top10 = beta + 1.645*sd	
gen sd_bot10 = beta - 1.645*sd
sum sd*
twoway rcap sd_top sd_bot age_centered, ///
	cmissing(n) lwidth(thin) lcolor(navy) lpattern(solid) xsc(r(-5 7)) xlabel(-5(1)7)  xline(-1) xline(3.5, lp(dash) ) ysc(r(-0.15 .15)) ylabel(-0.15(0.05).15,grid)   ///
	|| rcap sd_top10 sd_bot10 age_centered, cmissing(n) lwidth(thin) lcolor(navy) lpattern(solid) xsc(r(-5 7)) ///
	|| scatter beta age_centered, mcolor(navy) scheme(s1color) lcolor(navy) lpattern(solid) yline(0, lcolor(gs2)) ///
	title("Earnings, relative to age-49 average")  note("N=${n}") ///
	legend(off) ytitle("Estimated Coefficients") xtitle("Years Relative to Age 50")   
graph export "`dir_output'/figureD7c.png",  replace 	
restore

		
	
clear all
exit
 	

